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Abstract 

We propose a new analytical form of the quark-gluon plasma equation of state (EoS). This new EoS 
satisfies all qualitative features observed in the lattice QCD calculations and gives a good quantitative 
description of the lattice results in the SU(3) gluodynamics. The energy density for the suggested EoS 
looks similar to that in the bag model, but requires a negative value of the bag constant. 
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The transition from a confined hadron-resonance phase to a deconfined phase, the quark 
gluon plasma (QGP), is expected at high temperatures and/or baryonic chemical potentials. 
For several decades, the bag model (BM) equation of state (EoS) has been used to describe the 
QGP (see, e.g., Ref. l|). In the simplest form, i.e. for non- interacting massless constituents 
and zero values of all conserved charges, the BM EoS reads: 

e{T) = asBT' + B, p{T) = ^ T' - B , (1) 

where the energy density e and the pressure p have a simple dependence on temperature T 
modified by adding a positive constant B ("vacuum pressure"). The Stefan-Boltzmann (SB) 
constant asB in Eq. ([T]) equals to: 

^SB = ^ (dB + Idp^ , (2) 
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where (Ib and dp are the degeneracy factors for the bosons (gluons) and fermions (quarks and 
antiquarks), respectively. The zero value of the baryonic chemical potential in Eq. ([T]) is a valid 
approximation for the QGP created in nucleus-nucleus collisions at the BNL RHIC and even 
better for future experiments at the CERN LHC. Note that also most lattice QCD calculations 
for the QGP EoS correspond to zero or very small values of the baryonic chemical potential. 
Equation ([1]) is assumed to be valid at T > T^, where the critical temperature Tc corresponds 
to a 1*^* order phase transition in the pure SU(3) gluodynamics or to a smooth crossover in the 
full QCD. 

The main features of the QCD deconfined matter EoS can be illustrated by the Monte 
Carlo (MC) lattice results ^ for the SU(3) gluodynamics presented in Fig. [TJ They can be 
qualitatively summarized as follows: 

1. The pressure p{T) rapidly increases at T > T^, while at high T the system reaches the 
ideal massless gas behavior p = e/3. 

2. However, the constant a = e/T^ = Sp/T"^ observed at high T is about 10 20% smaller 
than the value of ctsb in Eq. ([2]). 

3. At high T, both e/T^ and 3p/T^ approach their limiting value a from below. 
Note that these properties of the gluon plasma EoS are also valid in the full QCD. 

The phenomenological success of the BM EoS ([1]) is due to the fact that it satisfies the first 
property: it gives p = e/3 at high T and shows an abrupt drop of p{T) near T^. However, the 
BM EoS is in a contradiction with the second and third features listed above (see Fig. [1]). The 
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FIG. 1: The MC lattice results for SU(3) gluodynamics [2] for the energy density (squares) and 
pressure (circles) at T > Tc- The size of the symbols corresponds to the error-bars reported in Ref. [^. 
The dashed horizontal line corresponds to the SB constant (JsB = Svr^/lS . The solid lines show the 
BM EoS ([1]) with d = 16 and B = 1.7T^ for e/T^ (upper line) and 3p/T'^ (lower line). 

BM EoS ([T]) gives no suppression of the SB constant. Note that we restrict our consideration 
to the present lattice results available at finite temperature interval < T < 4:.5Tc and do 
not discuss the possible asymptotic behavior at T — )• oo. The BM energy density e(T)/T^ 
approaches its SB limit from above. This contradicts the MC lattice results. Despite these 
evident problems, the BM EoS ([T]), due to its simplicity, is still one of the most popular models 
for phenomenological applications. 

In this letter we suggest a new analytical parametrization for the QGP EoS. It satisfies all 
three properties listed above, gives a good quantitative description of the MC lattice results for 
the SU(3) gluodynamics, and is almost as simple as Eq. ([T]). 

As the first step, we consider the suppression of the asB constant. For this purpose the 
quasi-particle approach of Ref. (see also recent papers and references therein) will be 
used. The system of interacting gluons is treated as a gas of non-interacting quasiparticles with 
gluon quantum numbers, but with mass m(T) which depends on T. The particle energy u and 
momentum k are assumed to be connected as u = [e + m2(T)]'/' . The energy density and 
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pressure take then the following form: 

= ^ / k'dk — - B\T) = po(T,a;) - B\T) , (4) 

where the degeneracy factor d = 2{N^ — 1) equals 16 for the SU(3) gluodynamics. The tem- 
perature dependent function B*{T) in Eq. ([3]) was introduced for the first time in Ref. It 
results from the thermodynamical relation, 

- PiT) = e{T) , (5) 



which leads to the equation for the function B*{T), 

dB* Ao(T,a;) dm 



(6) 
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where Aq = Eq — 3po, and Po defined by Eqs. fl3]H|) are the ideal gas expressions for 
massive bosons. If the function m{T) is known one can calculate B*{T) from Eq. ([6]) up to 
an arbitrary integration constant B. The linear relation m = aT with a = const > used 
for all T > Tc guarantees the high temperature behavior of e{T) and p{T) in agreement with 
the MC lattice results. For m = aT, the function B*{T) derived from Eq. equals to 
B*(T) = B — Ao(T, a;)/4 . One obtains the energy density ([3]) and the pressure (jlj), 

e{T) = aT' + B, p{T) = ^ T' - B , (7) 

where the modified SB constant a equals to: 
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The Ki and K2 in Eq. ([H]) are the modified Bessel functions. The constant a in Eq. includes 
the suppression factor n{a) which is defined by Eq. ((HI) and presented in Fig. [21 Therefore, 
an assumption of the linear T-dependent mass, m = aT, leads to the EoS similar to the 
bag model EoS, but with the suppressed SB constant ([8]). For a — > 0, k ^ 1 follows, and 
Eq. ([7]) coincides with Eq. ([1]). The modified SB constant a = 4.73 < asB allows to fit the high 
temperature behavior of s{T) and p{T). This requires /t(a) = 0.90 and a = 0.84. 

At the second step, which is the main point of our model construction, we include the linear 
in T contribution to the QGP pressure. If the function e{T) is known, Eq. ([5]) is a 1*^* order 




FIG. 2: The suppression factor n{a) from Eq. ([HI) as a function of the parameter o. 



differential equation for the function p{T). The general solution of this equation includes an 
arbitrary integration constant which results in the linear in temperature term in the function 
p{T). This was discussed for the first time in Ref. [s]. Thus, for e{T) in the form of Eq. (j7]), 
the general solution of Eq. (jSj) for p{T) can be written as follows, 



s{T) 



aT* + B 



p{T) 



1 
3 



B - AT . 



(9) 



A sum of the first and second terms in the expression for p{T) is a partial solution of the 
inhomogeneous differential equation ([5]) with s{T) given by ([9]), whereas the last term in p{T) 
corresponds to a general solution of the homogeneous equation Tdp/dT — p = 0. Therefore, the 
thermodynamical relation (|5]) between the pressure and energy density admits the linear in T 
contribution to p{T), which is fully invisible in the b{T) function. 

Equation ([9]) defines our model suggestion for the QGP EoS. For brevity we call it the A- 
bag model (A-BM). The formula for e{T) looks formally the same as in Eq. ([T]). However, 
the pressure function p{T) in the A-BM ([9]) contains one more parameter A comparing to 
the original BM EoS (II]). The model parameters, a = 4.73 and B = — 2.37 are found 
from fitting the MC lattice results {2] for the energy density function e{T). The third A- 
BM parameter A = 3.94 Tj} is fixed by fitting the pressure function p{T). One finds a good 
description of the MC lattice results for e{T) and p{T) within the A-BM EoS ([9]) for all ^ T > Tc 



^ To be precise, note that we consider the MC lattice points with T > 1.02 Tc to avoid the uncertainties of 
e(r) at T = Tc . 
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FIG. 3: The MC lattice results and the dotted horizontal line are the same as in Fig. 1. The solid lines 
correspond to the ^-BM EoS Q with a = 4.73, A = 3.94 T^, and S = - 2.37 for e/T'^ (upper 
line) and 3p/T^ (lower line). 

as shown in Fig. |3l The parameter a in Eq. (Q regulates the high temperature behavior of 
e/T^ ^ 3p/T^ ^ a. As A > 0, the hnear in T term gives a negative contribution to p{T) 
and guarantees both the correct high temperature asymptotic behavior of p{T) and its strong 
drop at T near T^. The bag parameter B in Eq. is found to be negative, in contrast to the 
positive bag constant B in the standard BM EoS ([1]). Thus, according to the A-BM (Q, e/T^ 
approaches its high temperature hmit a from below. This is in agreement with the MC lattice 
results. 

An important characteristic of the EoS is the so-called interaction measure, {e — 3p)/T^, 
which shows the deviation from the system of noninteracting massless particles. For the A-BM 
EoS dl]) the interaction measure reads, 

e — 3p 



rp4 



3A AB 



(10) 



The MC lattice results 



demonstrate a prominent maximum of the function [e — 3p)/T^ at 
Tmax — 1-1 Tc. The maximum of {e — 3p)/T^ is described in the A-BM iQ. This happens due 
to different signs of the A- and B-terms {A > 0, B < 0) in the r.h.s. of Eq. (ITOl) . Note that 
such a maximum is .produced ,y *e .o-called fu.y bag .ode. B. In .>7.„de., there 
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are contributions to both p{T) and e{T), 

e{T) = aT^ - CT^ + B , p{T) = - - C - B , (11) 

3 

and a comparison with the MC lattice results {2! gives C > and 5 > 0. In that case, 
[e — 2>p)/T'^ = 2C/T^ + AB /T^ corresponds to a monotonous decreasing function of T as both 
terms are positive. A comparison of the EoS ([9]) and f|TT]) will be discussed in more details 
in Ref. j^]- An extension of the A-BM to the SU(A^c) gluodynamics with A^^c > 3 [8], to the 
quark degrees of freedom and non-zero baryonic chemical potentials can be done along the same 
scheme and will be considered elsewhere. 

In summary, we have suggested a new EoS for the deconfined matter - the A-BM ([9]). It 
satisfies all qualitative features of the MC lattice results at T > Tc and gives a good quantitative 
description of the lattice results for the SU(3) gluodynamics, see Fig. [31 The expression for 
s{T) in the A-BM ^ looks similar to that in the BM ([T]). However, the pressure function 
p{T) in the A-BM ([9]) contains a new linear in T negative term which does not contribute to 
e{T). The presence of this negative pressure term leads to a principal difference between the 
bag term B in the BM and that in the A-BM. The bag parameter in the A-BM ([9]) is found 
to be negative, in contrast to the positive bag constant B in the BM EoS ([T]). The A-BM (j9]) 
gives a simple analytical parametrization of the QGP EoS. This opens new possibilities for its 
applications in the hydrodynamic description of the QGP. 
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